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ABSTRACT 

We consider MHD waves as a heating source of cool cores of galaxy clusters. 
In particular, we focus on transverse waves (Alfven waves), because they can 
propagate a longer distance than longitudinal waves (sound waves). Using MHD 
simulations, we found that the transverse waves can stably heat a cool core if the 
wave period is large enough (> 10 8 yr). Moreover, the longitudinal waves that 
are created as a by-product of the nonlinear evolution of the transverse waves 
could be observed as the 'ripples' found in cool cores. 

Subject headings: cooling flows — galaxies: active — galaxies: clusters: general 
— waves — X-rays: galaxies: clusters 



1. Introduction 

X-ray observations have shown that the radiative cooling time of the ICM (intracluster 
medium) in the central regions of galaxy clusters (cool cores) is generally smaller than the 
Hubble time. If there were no heating sources, the ICM flo ws subsonica lly toward the cluster 



center with a mass deposition rate of ~ 100-1000 M Q yr _1 (IFabianll 19941 ). This flow had been 
called a 'cooling flow'. However, recent X-ray observations have shown that the actual cooling 
rate of the ICM is much smaller. The Japanese ASCA (Advanced Satellite for Cosmology and 
Astrophysics) team indicated that metal emission lines from the low temperature cooling gas 
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were much weaker than had been predicted by the classical cooling flow mo del (llkebe et al. 



2001 



Kaastra et al 



2001 



1997; Makishima et al.l l200lf). This has been confirmed by XMM-Newton (IPeterson et al. 



Tamura et al.ll200ll ); the actual mass deposition rate is at least 5 



or 10 times less than that previously assumed. 

The lack of the cooling flows indicates that there are heating sources in the central 
regions of clusters that prevent the development of cooling flows. Most popular candidate 
of the heating source is AGNs at the centers of clusters. However, although the energy 
released by the AGNs is sufficient to cancel radiative cooling of the ICM in cool cores, the 
mechanism of energy transfer from the AGNs to the surrounding ICM is not clear. One 
possible mechanism is the dissipation of sound waves and we a k shocks that have evolved 
from sound waves (iFabian et al.l 120031 ; iRuszkowski et al.l 12004 ; iFabian et al.l 120051 ) . Those 
waves could be generated by the activities of the central AGNs of clusters. In fact, sound 
waves an d weak shocks are observed in the cen tral regions of the Perseus and the Virgo 
Clusters (iFabian et alfaoi librman et al.ll2005h . 



However, it has been shown that sound waves (longitudinal waves) and weak shocks 
that have evolved from the sound waves cannot stably heat a cool core, because they should 
dissipate before they propagate in the whole core. In other words, their dissipation length is 
too small. Even if their net energy is large enough to balance radiative cooling, the temper- 
ature profile of a cluster soon becomes irreg ular, which is not consistent with observations 
( IFujita &: Suzukill2005l ; iMathews et al.ll2006l ). On the other hand, it seems to be curious that 
structures that look like sound waves ('ripples') have been observe d far away (r ~ 50 kpc) 
from the central AGN in the Perseus cluster (IFabian et al.l 120061 ). Since the dissipation 
length of the waves is small, the density fluctuations or the ripples should not be observed 
there if they were created at the cluster center. 

So far, only longitudinal waves have been con sidered as heating sources in clusters . 
However, since the ICM is magnetized (~ 5-10 /zG; IClarke. Kronberg. &: Bohringerll200ll ). 
there should be magnetohydrodynamic (MHD) waves including tranverse waves. It has 
been known that tran sverse waves (Alfven waves) can propagate a longer distance than 



longitudinal waves (e. g ISuzukill2002l ) . Thus, the transverse waves could heat the whole core. 



In this letter, we consider the heating by tranverse waves and longitudinal waves that are 
a by-product of nonlinear evolution of the transverse waves. The longitudinal waves could 
be observed as the ripples. We do not restrict the wave generators only to AGNs, because 
a lo ng wave period is requir ed to reproduce the ripples as shown below. It is to be noted 
that iMakishima et al.l (120011 ) proposed a heating model of MHD effects, but they considered 
magnetic reconnection instead of waves. 
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Models 



The models presented here are based on heating models of the solar corona (jKudoh fc Shibata 



19991 ; ISuzuki fc Inutsukall2005l ). We study the propagation of waves along a magnetic field 
line extending from the cluster center in the radial direction. We assume axisymmetry about 
this field line. In the vicinity of the line, we consider a local orthogonal curvilinear coordinate 
system defined by r, the distance from the cluster center, by the azimuthal angle 9 measured 
about the axis of symmetry, and by £ measured in the direction perpendicular to both r and 
9 vectors (IHollweg. Jackson, fc Galloway! Il982l ). In the following, the indices, r, 9, and £, 
indicate the vector components for each direction. 

We assume ideal MHD (inviscid, perfect conductor) in a whole region: 
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where p is the ICM density, S is the cross section of the magnetic flux tube, v is the velocity, 
p is the pressure, B is the magnetic field, R is the distance from the axis of the symmetry, E 
is the electric field, j(— 5/3) is the adiabatic index, g is the gravitational acceleration, and 
L is the cooling rate. We assu med that 8/ 89 = 0, v% — 0, and B^ = in order to reduce 
the problem to one-dimension (IHollweg et al.lll982l ). For the sake of simplicity, we assume 
that the cluster is spherically symmetric, which means S oc r -2 . The initial magnetic field 
is given by B r oc r~ 2 and Bg = 0. We set the inner and outer boundaries at r in = 5 kpc 
and r out = 84 Mpc, respectively. Because of the large outer boundary, we do not have to 
consider the reflection of waves there. We adopted the relatively large inner radius to avoid 
the divergence of B r at the cluster center. The beta parameter of the plasma is set to be one 
at r = r; n . Thus, for r > r m it increases rapidly as r increases, which means that the Alfven 
velocity va is smaller than the sound velocity c s . At r = 20 kpc, we obtain B r = 6 p G. 
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Transverse waves are injected around the inner boundary in a form of acceleration in 
the 9 direction. The wave acceleration is given by 
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where gj n is the gravitational acceleration from the cluster potential at r 



and a and 



P are the parameters. Longitudinal waves are not injected for simplicity. It is to be noted 
that even if we give longitudinal waves with the same acceleration as that of the transverse 
waves, the resultant gas velocity of the longitudinal waves (v r ) is much smaller than that 
of the transverse waves (vg). This is because of the pressure gradient of the cluster, which 
prevents the ICM from being lifted in the r-direction. 



We adopt a cooling function based on the detailed calculations by ISutherland Dopita 



fll993h . 



L = n 2 A = [d{k B T) a + C 2 {k B Tf + C 3 ]mn 



(9) 



where n e and rij are the electron and the ion number densities, respectively, and the units 
for k B T are keV. For an average metallicity Z = 0.3 Z Q , the constants in equation (Q are 
a = -1.7, = 0.5, Ci = 8.6 x 10~ 3 , C 2 = 5.8 x 10~ 2 , and C 3 = 6.4 x 10" 2 , and we can 
approximate n,; n P = 0.70 (p/rriH) 2 , where ijih is the hydrogen mass. The units of A are 
10~ 22 ergs cm 3 (IRuszkowski fc Begelmanl 120021 1 . 



The model cluster we considered is similar to that adopted in § 3.2 of lFuiita. Suzuki, fc Wada 



(2004). The mass density profile of the cluster is given by the NFW profile (INavarro. Frenk. fc White 
19971 ). The viral mass is 1.2 x 10 15 M Q , the concentration parameter is 4.7, and the charac- 
teristic radius is 460 kpc. For this cluster, g m = 1300kms~ 1 (10 8 y r)~ 1 , which is a factor of two 

large r than the acceleration observed at the cold front in A1795 (jMarkevitch. Vikhlinin. &: Mazzotta 
20011 ). Initially, the cluster is isothermal and the temperature is 7 keV. The initial electron 
density at the cluster center is n e0 = 0.017 cm -3 . All of the calculations presented in this 
letter use CANS (Coordinated Astronomical Numerical Software )0. The number of grid 
points is 16392. The grid size is 8.3 pc for r < 120 kpc and it increases logarithmically for 
r > 120 kpc. 



3. Results 

First, we investigate the case where waves are not injected for comparison (a genuine 
cooling flow). Fig. [T] shows the evolution of temperature and electron number density as 
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functions of distance from the cluster center. The interval of the lines is 0.8 Gyr. Because of 
radiative cooling, the temperature decreases and the density increases monotonically at the 
cluster center. 

Next, we consider the heating by Alfven waves by injecting transverse waves at the 
inner boundary. We take a = 1 and P = 0.5 Gyr. For this wave period P, the wave 
length of the sound waves is c s P ~ 700 kpc. On the other hand, the wave length of the 
Alfven waves is smaller and is vaP ~ 10 kpc at r ~ 20-50 kpc, which is comparable to the 



observed wave length of ripples (IFabian et al.ll2006l ). This wave period may be too large to be 
attributed to the activities of the central AGN; it could be attributed to g as motion around 



the c entral galaxy i nduced by minor cluster mergers s uch as 'sloshing' (IMarkevitch et al. 



20011 ) or 'tsunami' (IFujita. Matsumoto. & Wadal 12004 ) . The wave injection in a form of 
acceleration (eq. [S]) would be preferable to model the waves of this kind, because the 
acceleration is independent of the gas flow toward the cluster center and the AGN. The 
transverse waves would be generated on the side of the sloshing central galaxy, but not in 
front and in the rear. 

Fig. [2] shows the evolution of the temperature and density profiles when transverse 
waves are considered. We found that while the temperature at the cluster center decreases 
for t < 4 Gyr, it is stable and almost constant for t > 4 Gyr until we stopped the calculation 
at t — 6 Gyr. Since we inject waves in a form of acceleration irrelevant of the gas mass of 
the core, the wave energy increases as the gas density in the core and thus the gas mass of 
the core increases. This is because for a given acceleration, the heavier the gas core is, the 
larger kinetic energy is given to the gas core. This is one reason that the temperature stops 
decreasing. The wave energy flux per unit area is given by Fa = pvj(vA + v r ) + v r Bj/(8n). 
At r = 15 kpc (outside the wave injection region), the energy flux is 4irr 2 FA ~ 10 42 erg s _1 
at t ~ 0, and it increases to ~ 10 44 erg s _1 at t ~ 6 Gy r. Numerical simulations s howed that 
the energy can be supplied by minor cluster mergers (ITittley fc Henriksenl 120051 ) . 



Fig. [3] shows the longitudinal and transverse velocity profiles at t — 1.6 and 3.2 Gyr. 
The transverse waves have reached r = 57 kpc at t — 1.6 Gyr and r = 66 kpc at t — 3.2 Gyr. 
In Fig. [3b, sharp jumps are seen; these jum ps are formed through the steeping of Alfven 
waves and are called 'Switch-on shocks' (e.g. lHoUweglll982l ; ISuzukill2004l ). At these shocks, 
while entropy is generated, longitudina l wave s are also generated through nonlinear co upling; 
(e.g. iHollweej Il992l ; iKudoh fc Shibatal Il999l : ISuzukil I2OO2I ; ISuzuki & Inutsukal l2005h . The 
longitudinal waves (sound waves) soon steepen and turn into weak shocks (N-waves; Fig. [3^), 
which also contribute to the heating of the ICM. Note that v r is much smaller than vg. 
Contrary to the transverse waves, the longitudinal waves can be observed in temperature 
and density profiles. The temperature and density fluctuations seen in Fig. [2] reflect the 
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generation of the longitudinal waves. In particular, the density fluctuations could be observed 
as ripples. Since these longitudinal waves are generated away from the cluster center, they 
can be observed even at r > 50 kpc (Fig. [2b). 

We also studied the mass inflow rate at the inner boundary. In the cooling flow case 
(Fig. HJ, the inflow rate at t = 4.8 Gyr is M = 220 M Q yr" 1 , and it is 310 M Q yr" 1 at 
t = 6 Gyr. On the other hand, waves can decrease the inflow rate significantly. In the case 
of Fig. [21 it is M = 60 M & yr _1 at t ~ 4.8 Gyr, which means that the waves effectively 
heat the surrounding ICM, although they cannot completely shut down a cooling flow. The 
inflow rate gradually increases, but even at t ~ 6 Gyr, it is still M = 80 M & yr _1 . 

We found that the heating by the Alfven waves is fairly stable even if we change the 
acceleration parameter a, because as we mentioned above, high density coming from radiative 
cooling results in a high heating rate, and vice versa. Fig. 0] shows the results when a = 2 
and P = 0.5 Gyr. Although the central temperature at which it stops decreasing (~ 4 keV; 
Fig. Ek) is higher than that in the case of a = 1 (~ 2 keV; Fig. [2^), the heating balances 
cooling after t ~ 4 Gyr. However, if a is much smaller than one, the waves cannot prevent 
the development of a cooling flow. In summary, a must be tuned to be ~ 1-2 in order to 
reproduce observed temperature profiles. 

Moreover, we have investigated the heating by Alfven waves with a relatively small wave 
period (P ~ 10 7 yr), which would be appropriate for waves generated by the activities of the 
central AGN. We found that the heating by these waves is not stable. The reason is similar 
to the heating by sound waves alone, that is, the waves can heat only the innermost region 
of the cluster. 



4. Summary and Discussion 

We have shown that Alfven waves (transverse waves) with a long wave period can 
effectively heat cool cores of galaxy clusters. Since the dissipation rate of the waves is small, 
the waves can reach the outer core and then heat the whole core before they dissipate. 
During the propagation, sound waves (longitudinal waves) are generated through nonlinear 
coupling, and they would be observed as ripples. The long dissipation length indicates the 
superiority of the heating by the Alfven waves over that by sound waves alone. In real 
clusters, however, it would be possible that sound waves created by AGN activities heat the 
inner cores, while Alfven waves induced by cluster mergers heat the outer cores. In this 
case, the wave acceleration (A w ) required to heat the cores will be much smaller than that 
we assumed. We assumed that the Alfven waves are generated close to the cluster center 
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(r ~ 5-10 kpc). iFujita et al.l (120041 ) showed that complex gas motion should be created 
even at the cluster center by minor cluster mergers. Thus, our assumption may be justified. 
However, if the waves are generated at some distance from the center, another heating source 
is required to heat the center. 

In this letter, we considered only spherically symmetric magnetic fields. If they are not 
spherically symmetric, the heating should spatially be biased. In the region with stronger 
magnetic fields, Alfven waves propagate faster and the heating is more effe ctive. This could 
result in multiphase gas observed in the Perseus cluster ( Fabian et al. 2006|). If the longevity 



of optical filaments observed in the Perseus cluster (IConselice. Gallagher, fc Wysd 120011 ) 
is smaller than the wave period, the filaments would indicate the direction of the wave 
oscillation or the motion of the central galaxy relative to the surrounding ICM. That is, 
the waves would be propagating in the direction perpendicular to the optical filaments. 
The model presented here predicts relatively large transverse velocity fluctuations (Fig. [3b), 
which could be detected in near-future space X-ray missions such as NeXT. 
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Fig. 2. — Same as Fig. [TJ but for a = 1 and P = 0.5 Gyr. 
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Fig. 3. — The longitudinal (v r ) and transverse (vg) velocity profiles at t — 1.6 and 3.2 Gyr. 
Note that the scales of the vertical axes are different. 




Fig. 4. — Same as Fig. [TJ but for a = 2 and P = 0.5 Gyr. 



